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The thermal conductivity of the color-flavor locked phase 
of dense quark matter is calculated. The dominant contri- 
bution to the conductivity comes from photons and Nambu- 
Goldstone bosons associated with breaking of baryon number 
which are trapped in the quark core. Because of their very 
large mean free path the conductivity is also very large. The 
cooling of the quark core arises mostly from the heat flux 
across the surface of direct contact with the nuclear matter. 
As the thermal conductivity of the neighboring layer is also 
high, the whole interior of the star should be nearly isother- 
mal. Our results imply that the cooling time of compact stars 
with color-flavor locked quark cores is similar to that of ordi- 
nary neutron stars. 

12.38.Aw, 21.65.+f 



I. INTRODUCTION 

At sufficiently high baryon density the nucleons in nu- 
clear matter should melt into quarks so that the system 
becomes a quark liquid. It should be weakly interacting 
due to asymptotic freedom jjj, however, it cannot be de- 
scribed as a simple Fermi liquid. This is due to the non- 
vanishing attractive interaction in the color antitriplet 
quark-quark channel, provided by one-gluon exchange, 
which renders a highly degenerate Fermi surface unsta- 
ble with respect to Cooper pairing. As a result the true 
ground state of dense quark matter is, in fact, a color 
superconductor 

Recent phenomenological |3| and microscopic studies 
[Q-Q have confirmed that quark matter at a sufficiently 
high density undergoes a phase transition into a color 
superconducting state. Phenomenological studies are ex- 
pected to be appropriate to intermediate baryon densi- 
ties, while microscopic approaches are strictly applicable 
at asymptotic densities where perturbation theory can 
be used. It is remarkable that both approaches concur 
that the superconducting order parameter (which deter- 
mines the gap A in the quark spectrum) lies between 10 
and 100 MeV for baryon densities existing in the cores 
of compact stars. For a review of color superconducting 
quark matter see Ref . || . 

At realistic baryon densities only the three lightest 
quarks can participate in the pairing dynamics. The fact 
that the current mass of the strange quark is more than 
an order of magnitude larger than the masses of the up 
and down quarks suggests two limiting cases. The first 



extreme limit assumes that the strange quark is suffi- 
ciently massive that it does not participate in the pair- 
ing dynamics at all. The corresponding quark phase is 
the color superconducting phase involving two quark fla- 
vors (2SC). Since it was suggested in Ref. || that such a 
phase is absent in compact stars, we focus on the second 
extreme case. This applies when all three quark flavors 
have a mass which is small in comparison to the chemi- 
cal potential so that they participate equally in the color 
condensation. The ground state is then the so-called 
color flavor locked (CFL) phase jn|. Here the original 
gauge symmetry SU(3) C and the global chiral symme- 
try SU(3)l x SU(3)r break down to a global diagonal 
"locked" SU(3) c+ l+r subgroup. Because of the Higgs 
mechanism the gluons become massive and decouple from 
the infrared dynamics. The quarks also decouple because 
large gaps develop in their energy spectra. The break- 
ing of the chiral symmetry leads to the appearance of an 
octet of pseudo-Nambu-Goldstone (NG) bosons (ir° , tt , 
K^, K°, K°, rj). In addition an extra NG boson </> and a 
pseudo-NG boson 77' appear in the low energy spectrum 
as a result of the breaking of global baryon number sym- 
metry and approximate U(1)a symmetry, respectively. 

The low energy action for the NG bosons in the 
limit of asymptotically large densities was derived in 
Refs. HQ. By making use of an auxiliary "gauge" 
symmetry, it was suggested in Ref. [[l3| that the low en- 
ergy action of Refs. §T],[l^] should be modified by adding 
a time-like covariant derivative to the action of the com- 
posite field. Under a favorable choice of parameters, the 
modified action predicted kaon condensation in the CFL 
phase. Some unusual properties of such a condensate 
were discussed in Refs. [|14|]15| 1. 

While the general structure of the low energy action 
in the CFL phase can be established by symmetry argu- 
ments alone ]rjj] , the values of the parameters in such an 
action can be rigorously derived only at asymptotically 
large baryon densities p^ , ^3[ . Thus, in the most inter- 
esting case of intermediate densities existing in the cores 
of compact stars, the details of the action are not well 
known. For the purposes of the present paper, however, 
it suffices to know that there are 9 massive pseudo-NG 
bosons and one massless NG boson cf> in the low energy 
spectrum. If kaons condense Jlj| an additional NG boson 
should appear .[| These NG bosons should be relevant for 



1 Strickly speaking the additional NG boson in a phase with 
kaon condensation has a small, but nonzero, mass as a result 
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the kinetic properties of dense quark matter. 

It has been found |l(if| that neutrino and photon emis- 
sion rates for the CFL phase are very small so that they 
would be inefficient in cooling the core of a neutron star. 
The purpose of the present investigation is to determine 
quantitatively the thermal conductivity of the CFL phase 
of dense quark matter in order to see whether it can sig- 
nificantly impact the cooling rate. We shall argue that 
the temperature of the CFL core, as well as the neighbor- 
ing neutron layer which is in contact with the core, falls 
quickly due to the very high thermal conductivities on 
both sides of the interface. In fact, to a good approxima- 
tion, the interior of the star is isothermal. A noticeable 
gradient of the temperature appears only in a relatively 
thin surface layer of the star where a finite flux of energy 
is carried outwards by photon diffusion. 

The paper is organized as follows. In the next section, 
we derive the microscopic expression for thermal conduc- 
tivity of the (pseudo-) NG bosons in the CFL phase. In 
Sec. [II we estimate the mean free path £ of the NG boson 
in the CFL phase. Knowledge of this quantity is crucial 
to a reliable calculation of the thermal conductivity. The 
role of photons in affecting the thermal conductivity and 
specific heat of the quark core is discussed in Sec 
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VL In the Appendix we present some useful for- 



In Sec. [v| we discuss our findings and their implication 
for the cooling dynamics of a compact star with a quark 
core in the CFL phase. Our conclusions are summarized 
in Sec. 
mulae. 



II. THERMAL CONDUCTIVITY 

A detailed understanding of the cooling mechanism of 
a compact star with a quark core is not complete with- 
out a study of thermal conductivity effects in the color 
superconducting quark core. The conductivity, as well as 
the other kinetic properties of quark matter in the CFL 
phase, is dominated by the low energy degrees of free- 
dom. It is clear then that at all temperatures of interest 
to us T <C A, it is crucial to consider the contributions 
of the NG bosons. In addition, there may be an equally 
important contribution due to photons; this is discussed 
in Sec. 



IV, Note that, at such small temperatures, the 



gluon and the quark quasiparticles become completely ir- 
relevant. For example, a typical quark contribution to a 
transport coefficient would be exponentially suppressed 
by the factor exp(— A/T). Therefore, in the rest of this 
section, we concentrate exclusively on the contributions 
of the NG bosons. 

Let us start from the general definition of the ther- 
mal conductivity as a characteristic of a system which is 



of weak interactions that break strangeness. The estimated 
value of the mass is of order 50 keV fll5|]. 



forced out of equilibrium by a temperature gradient. In 
response to such a gradient transport of heat is induced. 
Formally this is described by the following relation: 



= -ndjT, 



(1) 



where Ui is the heat current, and k is the heat conduc- 
tivity. As is clear from this relation, the heat flow would 
persist until a state of uniform temperature is reached. 
The higher the conductivity, the shorter the time for this 
relaxation. 

In the linear response approximation, the thermal con- 
ductivity is given in terms of the heat current correlator 
by a Kubo-type formula. We derive the expression for 
the heat (energy) current carried by a single (pseudo-) 
NG boson field tp. The corresponding Lagrangian den- 
sity reads 



L = | (d Q (pd (p - v 2 di(pdiip - m 2 ip 2 ) + 



(2) 



where the ellipsis stand for the self-interaction terms as 
well as interactions with other fields. Notice that we 
introduced explicitly the velocity parameter v. In mi- 
croscopic studies of color superconducting phases, which 
are valid at very large densities, this velocity is equal to 
l/\/3 for all (pseudo-) NG bosons. It is smaller than 
I because Lorentz symmetry is broken due to the finite 
value of the quark chemical potential. By making use 
of the above Lagrangian density, we derive the following 
expression for the heat current: 



i)L 



d ip = v 2 diipd tp, 



(3) 



This definition leads to the expression |17j for the heat 
conductivity in terms of the corresponding correlator: 

where, in the Matsubara formalism, 

I!jj(zf2 m ) = v T j ^ J -j^—^kikjiVl n {iQ, n +if2 n _ m ) 

x S(iO, n , k)S(iQ n - m , k). (5) 

Here Q n = 2-ktiT is the bosonic Matsubara frequency, 
and S(iQ n , k) is the propagator of the (pseudo-) NG bo- 
son. In general, the propagator should have the following 
form: 



S(u,k) = 



1 



(cu + iT/2) 2 -v 2 k 2 - 



(6) 



where the width parameter r(w, k) is related to the in- 
verse lifetime (as well as the mean free path) of the boson. 

In our calculation, it is very convenient to utilize the 
spectral representation of the propagator, 
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S(iQ n , k) = - 



du>A(u>, k) 
zf2„ — LU 



(7) 



Then, the conductivity is expressed through the spectral 
function A(u>,k) as follows: 



2ttT 2 



uj 2 duj 



sinh z 



2T 



d 3 k 



hk 3 A 2 (uj,k). (8) 



By making use of the explicit form of the propagator in 
Eq. (^), we see that the spectral function of the (pseudo-) 
NG boson is 



A(u,k) 



CUT 



(cj 2 - e 2-r 2 /4) 2 + w 2 r 2 



r 

4e fc 



{LU-e k ) z +T 2 /A 



— (lu —>■ —lu) 



(9) 



where e k = \Jv 2 k 2 + m 2 . By substituting this spectral 
function into the general formula (|J), we notice that the 
off diagonal components of the thermal conductivity are 
zero and all the diagonal components are equal. Both of 
these facts are a consequence of the rotational symme- 
try of the system. Thus the conductivity is characterized 
by a single scalar quantity k which is introduced as fol- 
lows: Kij = n5ij . The explicit expression for this scalar 
function reads 

1 f°° ujduj 



K 



48\/27r 2 wrT 2 J sinh 



°— (Vx 2 + iu 2 r 2 + x 

2T 



3/2 



(10) 



where we introduced the notation X = lu 2 — m 2 — T 2 /4. 
In the calculation leading to Eq. ([l0|) , we made use of the 
result p): 



c^x^dx 



[(a; 2 + 6) 2 + a 2 ^ g 



V2|c 



2 + b 2 -b 



3/2 



(11) 



For our purposes it will be sufficient to consider the con- 
ductivity in the limit of a very small width. This is be- 
cause the (pseudo-) NG bosons in the CFL quark matter 
are weakly interacting. In this case we derive the follow- 
ing approximate relation: 



1 



24n 2 vT 2 T 



(Lulu 



sinh^ 



(- 2 



2\3/2 



(12) 



2T 



It is interesting to notice that the same result can be 
easily derived from Eq. (@). Indeed, in the limit of a very 
small width, i.e., r — * 0, we can replace the square of the 
spectral function by a sum of two (5-functions: 



A 2 (to,k) 



4r4 



[<5 (u - e k ) + 5(w + e k )] 



(13) 



where the coefficients in front of the (5-functions were de- 
termined unambiguously by using the following replace- 
ment rule: 



lim -r^, „ ; . „ — > AttS(x) 

r^o (x 2 + T 2 /4) 2 v ' 



(14) 



Now in the limit of small temperature, T <^ m, the result 
in Eq. ( [l2|) is given approximately by: 



m 



x dx 



2Att 2 vTT 2 Ju 

m 5 / 2 VT 
2 % /2 7 r 3 / 2 wr f 



sinh 



m/T 



2 7re\/l+ic 2 



IT 



(15) 



This demonstrates clearly that the contributions of 
heavy pseudo-NG bosons to the thermal conductivity 
are strongly suppressed. The largest contribution comes 
from the massless NG boson <f> for which the thermal 
conductivity is 



4T 3 



37r 2 wr 



x dx 



2tt 2 T 3 



<t> Jo 



sinh a; 45tT^ 



(16) 



In order to calculate the width is required, or 
equivalently the mean free path since 1$ = v/Y$, 
where v is the average thermal velocity of the particles 
responsible for heat transfer. This will be discussed in 
the next section. For the moment we note that expres- 
sion ( |T6| ) for the thermal conductivity is similar to the 
phonon conductivity often used in solid state physics [fl9[ . 
It is also interesting to note that the result for the ther- 
mal conductivity of both massive and massless bosons in 
Eqs. ( |l5| ) and ([If]) are consistent with a general relation 
from classical Boltzmann kinetic theory fl2Cfl: 



(17) 



where c v is the specific heat. In order to see that this 
relation holds, we note that the average thermal velocity, 
(ctefc / dk) , of massive bosons is a factor y/3T/m smaller 
than the velocity of massless bosons v — l/v3. After tak- 
ing this into account, as wel l as us ing t he expressions for 
the specific heat in Eqs. (A.8) and ( A12j ) in the Appendix, 
we find that the relation (p_7[) between the conductivity 
and the specific heat is indeed satisfied. 



III. MEAN FREE PATH OF THE NG BOSON 

In the preceding section we derived a general expres- 
sion for the boson portion of the thermal conductivity. 
The result involved the decay width of the (pseudo-) NG 
boson or equivalently the mean free path, £ = v/T, for 
which we derive a simple estimate here. Using this we 
can understand the effect of the heat conductivity on the 
relaxation time of the temperature gradient inside the 
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quark core of a star. This is crucial for understanding the 
cooling mechanism of compact stars with quark cores. 

As we have remarked, the contribution of massive 
pseudo-NG bosons to the thermal conductivity is sup- 
pressed by the exponential factor exp(— m/T). In the 
CFL phase of quark matter, however, there is one truly 
massless NG boson which should therefore give the 
dominant contribution to the heat conductivity. The in- 
teractions of <fi with the CFL matter leads to a finite value 
for its mean free path. Since this boson is a composite 
particle there is always a non-zero probability at finite 
temperature for its decay into a pair of quark quasiparti- 
cles. It is natural to expect that such a process is strongly 
suppressed at small temperatures, T« A. This is con- 
firmed by a direct microscopic calculation in the region 
of asymptotic densities which yields a decay width [Ell: 



"/</ 



(k) 



4(21 -81n2) 



exp 




(18) 



If this were the only contribution, then the order of mag- 
nitude of the mean free path of the NG boson would be 



T 



exp 




(19) 



This grows exponentially with decreasing the tempera- 
ture. For example, if T < A/33 the mean free path is 
30 km or more (in deriving this estimate we set A ~ 50 
MeV). This scale is a few times larger than the typical 
size of a compact star. The underlying physics in this 
calculation of the mean free path closely resembles the 
propagation of sound waves in superfluid helium. We re- 
call that there are two types of excitations in superfluid 
helium: gapless phonons and finite gap rotons. The main 
contribution to the thermal conductivity comes from the 
phonons whose mean free path is mainly determined by 
scattering on rotons ]22| . The value of the mean free 
path, as in the case of NG bosons, becomes exponentially 
large at small temperatures (T <C A). 

Now, in the case of the NG bosons the decay channel 
into quarks is not the only contribution to the mean free 
path since the NG bosons can also scatter on one another. 
The amplitude has been derived by Son |23| and is of 
order fc 4 //x 4 which gives a cross section of ~ T 6 //j, 8 . 
This yields the following contribution to the width of the 
NG boson: 



r rf 



ji9 



(20) 



where is the equilibrium number density of the NG 
boso ns. The explicit expression for is given in 
Eq. (All) in the Appendix. At small temperatures the 



scattering contribution in Eq. ( pp| ) dominates the width. 
This leads to a mean free path 



£r ps 8 x 10 ; 



^500 
± MeV 



km. 



(21) 



Here we defined the following dimensionless quantities: 
^soo = M /(500 MeV), A 50 = A/(50 MeV), and T MeV = 
T/(l MeV). Both and l^ qq depend very strongly 
on temperature, however the salient point is that they 
are both larger than the size of a compact star for tem- 
peratures TMeV of order 1. 

We define T to be the temperature at which the mas- 
sive NG bosons decouple from the system. This is deter- 
mined by the mass of the lightest pseudo-NG boson for 
which it is not presently possible to give a reliable value. 
Different model calculations |l2||l3|,Q produce different 
values which can range as low as 10 MeV. Thus, con- 
servatively, we choose T ~ 1 MeV. Then the mean free 
path of the NG boson is comparable to or even larger 
than the size of a star for essentially all temperatures 
T < T. It is also important to note that the mean free 
path is very sensitive to temperature changes. In par- 
ticular, at temperatures just a few times higher than T 
the value of I may already become much smaller than 
the star size. This suggests that, during the first few 
seconds after the supernova explosion when the temper- 
atures remain considerably higher than T, a noticeable 
temperature gradient may exist in the quark core. This 
should relax very quickly because of the combined effect 
of cooling (which is very efficient at T ^> T) and dif- 
fusion. After that almost the whole interior of the star 
would become isothermal. 

Before concluding this section, we point out that the 
geometrical size of the quark phase limits the mean free 
path of the NG boson since the scattering with the 
boundary should also be taken into account. It is clear 
from simple geometry that I ~ i?0; where Rq is the ra- 
dius of the quark core. A similar situation is known to 
appear in high quality crystals at very low temperatures 



IV. PHOTON CONTRIBUTIONS 

In this section, we discuss the role of photons in the 
CFL quark core. It was argued in Ref. |l(| that the mean 
free path of photons is larger than the typical size of a 
compact star at all temperatures T < T. One might 
conclude therefore that the photons would leave a finite 
region of the core in a very short amount of time after the 
core becomes transparent. If this were so, the photons 
would be able to contribute neither to thermodynamic 
nor to kinetic properties of the quark core. However the 
neighboring neutron matter has very good metallic prop- 
erties due to the presence of a considerable number of 
electrons. As is known from plasma physics, low fre- 
quency electromagnetic waves cannot propagate inside a 
plasma. Moreover, an incoming electromagnetic wave is 
reflected from the surface of such a plasma [2f| . In par- 
ticular, if fl p is the value of the plasma frequency of the 
nuclear matter, then all photons with frequencies lu < f2 p 
are reflected from the boundary. This effect is similar to 
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the well known reflection of radio waves from the Earth's 
ionosphere. 

The plasma frequency is known to be proportional to 
the square root of the density of charge carriers and in- 
versely proportional to the square root of their mass. 
It is clear therefore that the electrons, rather than the 
more massive protons, will lead to the largest value of 
the plasma frequency in nuclear matter. Our estimate 
for the value of this frequency is 



m e m p 



4.7 x 10 2 W — MeV, (22) 
V Po 



where the electron density n e = Y e p/m p is given in terms 
of the nuclear matter density p and the proton mass m p . 
Also m e denotes the electron mass, Y e ~ 0.1 is the num- 
ber of electrons per baryon, and po ~ 2.8 x 10 14 g cm" 3 
is equilibrium nuclear matter density. 

Since £l p is more than 100 MeV electromagnetic waves 
of essentially all frequencies for temperatures appropriate 
to the CFL phase will be reflected back into the core 
region. Thus photons present very early in the life of the 
star, those produced by decays in the CFL matter as well 
as those produced at the nuclear interface will be trapped 
in the core. In a way the boundary of the core looks like 
a good quality mirror with some leakage which will allow 
a thermal photon distribution to build up, even at the 
relatively high temperatures that existed during the first 
moments of stellar evolution. 

Now, since photons are massless they also give a siz- 
able contribution to the thermal conductivity of the CFL 
phase. The corresponding contribution k 7 will be similar 



to the contribution of the massless NG boson in Eq. (16). 
Since the photons move at approximately the speed of 
light p6| (v ~ 1 in our notation) and they have two po- 
larization states, we obtain 



4tt 2 T 3 
45I\ ' 



V. STELLAR COOLING 

In discussing the cooling mechanism for a compact star 
we have to make some general assumptions about the 
structure of the star. We accept without proof that a 
quark core exists at the center of the star. The radius 
of such a core is denoted by Rq, while the radius of the 
whole star is denoted by R. We exclude the possibility 
that the star is made completely of CFL quark matter. 
Schematically, the internal structure of the star is shown 
m Fig. 0. The quark core stays in direct contact with the 
neighboring nuclear matter. From this nuclear contact 
layer outwards the structure of the star should essentially 
be the same as in an ordinary neutron star. 




(23) 



FIG. 1. The schematic structure of a compact star. At 
the center is a quark core which is surrounded by ordinary 
nuclear matter. On the outside is a crust of nuclei, neutrons 
and electrons. 



Since the thermal conductivity is additive the total 
conductivity of dense quark matter in the CFL phase 
is given by the sum of the two contributions: 



2tt 2 , 

KCFL = + K 7 ~ — — T i? , 



(24) 



where for both a photon and a NG boson the mean free 
path I ~ Rq. This yields the value 

kcfl ^ 1-2 x 10 32 T 3 /ey i? , fcm erg cm-W^KT 1 , (25) 

which is many orders of magnitude larger than the ther- 
mal conductivity of regular nuclear matter in a neutron 
star Wh. 



A detailed analysis of the interface between the quark 
core and the nuclear matter was made in Ref . ]2§| ] . A sim- 
ilar analysis might also be very useful for understanding 
the mechanism of heat transfer from one phase to the 
other here. For our purposes it is sufficient to know only 
that direct contact between the quark and nuclear phases 
is possible. In fact the only assumption we need is that 
the temperature is continuous across the interface. 

Having spelled out the basic assumptions, we can con- 
sider the physics that governs stellar cooling. Let us start 
from the moment when the star is formed in a supernova 
explosion. Right after the explosion a lot of high-energy 
neutrinos are trapped inside the star. After about 10 to 
15 seconds most of them escape from the star by diffu- 
sion. The presence of the CFL quark core could slightly 
modify the rate of such diffusion ]29|-|3l|] . By the end of 
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the deleptonization process, the temperature of the star 
would rise to a few tens of MeV. Then, the star relatively 
quickly cools down to about T by the efficient process 
of neutrino emission. It is unlikely that the quark core 
would greatly affect the time scale for this initial cooling 
stage. An ordinary neutron star would continue to cool 
mostly by neutrino emission for quite a long time even 
after that |52| . Here we need to discuss how the presence 
of the CFL quark core affects the cooling process of the 
star after the temperature drops below T, 

Our result for the mean free path of the NG boson 
demonstrates that the heat conductivity of dense quark 
matter in the CFL phase is very high. For example, a 
temperature gradient of 1 MeV across a core of 1 km in 
size is washed away by heat conduction in a very short 
time interval of order km /v£(T) ~ 6 x 10~ 4 sec. In 
deriving this estimate, we took into account the fact that 
the specific heat and the heat conductivity in the CFL 
phase are dominated by photons and massless NG bosons 
and that £ ~ Rq. In addition, we used the relation in 
Eq. (|l7|). This estimate proves that, to a good approxi- 
mation, the quark core is isothermal at all temperatures 
T <f. 

The heat conductivity of the neighboring nuclear mat- 
ter is also known to be very high because of the large 
contribution from degenerate electrons which have a very 
long mean free path. It is clear, then, that both the 
quark and the nuclear layers should be nearly isothermal 
with equal values of the temperature. When one of the 
layers cools down by any mechanism, the temperature 
of the other will adjust almost immediately due to the 
very efficient heat transfer on both sides of the interface. 
This observation has a very important consequence for 
the cooling rate of the star. It will eventually be de- 
termined by the combination of two effects: (i) neutrino 
emission from the nuclear part of the star, and (ii) dif- 
fusion of photons through the outer crust made of nuclei 
and non-degenerate electrons. Even if, for any reason 
the intermediate nuclear layer happens to be abnormally 
thin photon diffusion will provide efficient cooling. The 
corresponding cooling time would be of order 10 7 years. 
This is many orders of magnitude shorter than the es- 
timated time of cooling of the CFL matter by neutrino 
emission Jl6[ . 

Now consider the order of magnitude of the cooling 
time for a star with a CFL quark core. One of the most 
important components of the calculation of the cooling 
time is the thermal energy of the star which is the amount 
of energy that is lost in cooling. There are contributions 
to the total thermal energy from both the quark and the 
nuclear parts of the star. The dominant amount of ther- 
mal energy in the CFL quark matter is st ored in photons 
and massless NG bosons, see Eq. ( AI2| ) and the brief 
discussion in the Appendix. It is given by the following 
expression: 



E, 



< ,?L(T) = ^-£c v (T')dT> 



6(l + 2v 3 )T fnTRo^ 



3v 



(26) 



where, for simplicity, we use the Newtonian approxima- 
tion. The specific heat here contains both <p and photon 



contributions 



which can be obtinaed 



from Eq. (A12) using the appropriate degeneracies and 
velocities. Numerically the expression in Eq. (^6|) leads 
to the following result: 



Ecfl{T) 



2.1 x 10 42 i?; 



O.fcm^MeV er g- 



(27) 



Here Ro,km is the quark core radius measured in kilo- 
meters. The thermal energy of the outer nuclear layer 
is provided mostly by degenerate neutrons. A numerical 
estimate is S: 



E NM (T) ~ 8.1 X 10 



1!) 



M - Mi 



Mr: 



^) 2/3 T\ leV erg, (28) 



where M is the mass of the star, Mo is the mass of the 
quark core and M Q is the mass of the Sun. It is cru- 
cial to note that the thermal energy of the quark core 
is negligible in comparison to that of the nuclear layer. 
Moreover, as the star cools the ratio Ecfl/Enm wm 
further decrease. 

The second important component that determines stel- 
lar cooling is the neutrino and/or photon luminosity 
which describes the rate of energy loss. Typically, the 
neutrino luminosity dominates the cooling of young stars 
when the temperatures are still higher than about 10 keV 
and after that the photon diffusion mechanism starts to 
dominate. As argued in Ref. p6| , neutrino emission from 
the CFL quark phase is strongly suppressed at low tem- 
peratures, i.e., at T < T in our notation. The neigh- 
boring nuclear layer, on the other hand, emits neutrinos 
quite efficiently. As a result, it cools relatively fast in the 
same way as an ordinary neutron star. The nuclear layer 
should be able to emit not only its own thermal energy, 
but also that of the quark core which constantly arrives 
by the very efficient heat transfer process. The analysis 
of this cooling mechanism, however, is greatly simplified 
by the fact that the thermal energy of the quark core is 
negligible compared to the energy stored in the nuclear 
matter. By making use of the natural assumption that 
local neutrino emissivities from the nuclear matter are 
not affected by the presence of the quark core, we con- 
clude that the cooling time of a star with a quark core 
is essentially the same as for an ordinary neutron star 
provided that the nuclear layer is not extremely thin. 



VI. CONCLUSIONS 

As is clear from our analysis, the thermal conductivity 
of the CFL color superconducting dense quark matter is 
very high for typical values of the temperature found in 
a newborn compact star. This is a direct consequence of 



G 



the existence of the photon and the massless NG boson 
associated with the breaking of baryon number. (Note 
that the quark contribution to the thermal conductivity 
is strongly suppressed in the CFL phase, as discussed in 
Rcf. [p4| .) The mean free path of the NG boson appears 
to be very large as a result of its weak self-interaction 
and its weak interaction with the rest of the matter in 
the CFL phase. The same is true for the photon. Our 
rough estimate shows that even at a relatively early stage 
of stellar cooling when the temperature is of order T, say 
1 MeV or so, the mean free paths are already larger than 
the size of the star. It should be noted, however, that 
the NG bosons cannot leave the CFL phase and escape 
from the core of the star carrying away their energy. By 
their nature, these NG bosons are collective excitations 
of the CFL phase and thus they are "confined" to the 
medium which made their existence possible. In this re- 
spect, they resemble phonons in solid state physics which 
also, by their nature, cannot escape from inside a crys- 
tal. Although for a different reason, thermal photons also 
cannot leave the finite region of the quark core. They are 
efficiently reflected from the electron plasma of the neigh- 
boring nuclear matter. 

We mention that the (pseudo-) NG bosons and pho- 
tons should also dominate other kinetic properties of 
dense quark matter in the CFL phase. For example, the 
shear viscosity should be mostly due to photons and the 
same massless NG bosons associated with the breaking 
of baryon number. The electrical conductivity, on the 
other hand, would be mostly due to the lightest charged 
pseudo-NG boson, i.e., the K + . Thus, in the limit of 
small temperatures T — > the electrical conductivity will 
be suppressed by a factor exp(— m K +/T). 

Since the neutrino emissivity of the CFL core is 
strongly suppressed, the heat is transferred to the outer 
nuclear layer only through the direct surface contact. 
While both the core and the outer layer contribute to the 
heat capacity of the star, it is only the outer layer which 
is capable of emitting this heat energy efficiently in the 
form of neutrinos. The combination of these two factors 
tends to extend the cooling time of a star. However, be- 
cause of the very small thermal energy of the quark core 
compared to that of the nuclear matter, the time scale 
for cooling would only get noticeably shorter than that 
for an ordinary neutron star when the quark core radius 
is nearly the same as the stellar radius. Thus it appears 
that the cooling of stars with CFL quark cores will differ 
little from the cooling of standard neutron stars. A sim- 
ilar conclusion has been reached for stars with regular, 
non-CFL quark interiors fl35| . 

In passing it is interesting to speculate about the possi- 
bility that a bare CFL quark star made entirely of dense 
quark matter could exist. If it were possible, it would 
look like a transparent dielectric |36f| . Our present study 
suggests such a star would also have very unusual thermal 
properties. Indeed, if the star has a finite temperature 
T < T after it was created, almost all of its thermal en- 
ergy would be stored in the NG bosons. Notice that all 



the photons would leave the star very soon after trans- 
parency set in because the star is assumed to have no 
nuclear matter layer. The local interaction as well as the 
self-interaction of the NG bosons is very weak so that 
we argued in Sec. Ill that their mean free path would 
be limited only by the geometrical size of the star. This 
suggests that the only potential source of energy loss in 
the bare CFL star would be the interaction of the NG 
bosons at the star boundary, together with very strongly 
suppressed jl(| photon and neutrino emission. It is likely, 
therefore, that such stars might be very dim and might 
even be good candidates for some of the dark matter in 
the Universe. A more detailed discussion of this issue is, 
however, beyond the scope of the present paper. 

Finally, we would also like to note that the thermal 
conductivity of the 2SC color superconducting phase of 
two lightest quarks should also be very large, but for 
a different reason. In this case there are quarks which 
do not participate in color condensation and thus give 
rise to gapless quasiparticles. The corresponding contri- 
bution to the heat conductivity (which is similar to the 
contribution of degenerate electrons inside nuclear mat- 
ter) should be very large. Therefore, if it were possible 
for the core of a compact star to be made of the 2SC 
quark phase, it would again be nearly isothermal. The 
results of Ref. || suggest, however, that the 2SC phase 
cannot appear in a compact star. 
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APPENDIX A: THERMODYNAMIC AL 
QUANTITIES 

In this Appendix, for convenience we list expressions 
for the thermodynamic quantities that we use throughout 
the main text of the paper. 

It is most convenient to start with the expression for 
the pressure. The other quantities (entropy, free energy, 
etc.) can then be easily derived In the case of a 

massive boson field, the pressure is 



p(bos) 
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d 3 k 
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-e k /T 
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where = \/v 2 k 2 + to 2 . The corresponding energy den- 
sity, number density and the specific heat are: 
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In the case of a small temperature (as compared to the 
value of mass), these thermodynamic quantities approach 
the following asymptotic expressions |pq| : 
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And, in the limit of massless particles, they are 
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where the value of the Ricmann function £(3) » 1.202. 
For completeness, let us also present some thermodynam- 
ical properties of color superconducting quarks. In this 
case, it is very convenient to separate the zero and finite 
temperature contributions explicitly. For the pressure of 
a single quark with a gap A, we find: 
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where E^' — y/ (e^ ± /i) 2 + A 2 are the quasiparticle en- 
ergies and /i is the chemical potential. The zero temper- 
ature contribution is approximately given by 
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The last term is needed so that the gap equation is ob- 
tained when the pressure is maximized with respect to 
variations in A. Using the gap equation to eliminate a 
and renormalizing by removing the divergent term, the 
pressure becomes 



(0) 



qrk 



127T 2 



M 2 A 2 
4tt 2 



in agreement with the results obtained using the frame- 
work of the Nambu-Jona-Lasinio model |3(| and directly 
from quantum chromodynamics p9[ . The ellipsis de- 
notes sublcading terms suppressed by powers of cither 
(m/» 2 or (A/V) 2 . As is easy to check, at small temper- 
atur es, i.e ., T -c A, the finite temperature corrections in 
Eq. ( |A13| ) are suppressed by a factor exp(— A/T). Thus, 
the above zero temperature contribution will dominate 
the pressure of the CFL phase. Similarly, the main con- 
tribution to the number density and the energy density 
of CFL matter at T -C A com es from the zero tempera- 
ture terms. Using Eq. ( A15Q the T = contribution to 
the number density is 



l qrk 



3tt 2 



1 + 



£ OA 
A 9^7 
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Note that the partial derivative term is needed since the 
gap equation, which specifies A(/x), was used in arriving 
at Eq. JA15[ ). This term contributes at the same or- 
der (unity) as the first term in parentheses. The energy 
density is then 



U, 



(0) 



qrk 



4tt 2 



M 2 A 2 
4tt 2 
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2/x dA 
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Here it is appropriate to note that all thermodynamic 
functions which are defined through the temperature 
derivatives of either the pressure (e.g., entropy) or the 
energy (e.g., specific heat) are exponentially suppressed 
at T <C A. Therefore, the dominant contributions to this 
second type of functions would come from (pseudo-) NG 
bosons in the CFL phase of dense quark matter. 
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